GENERALIZED SOLUTIONS OF THE VLASOV-POISSON 
SYSTEM WITH SINGULAR DATA 
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Abstract. We study spherically symmetric solutions of the Vlasov-Poisson 
system in the context of algebras of generalized functions. This allows to model 
highly concentrated initial configurations and provides a consistent setting for 
studying singular limits of the system. The proof of unique solvability in our 
approach depends on new stability properties of the system with respect to 
perturbations. 



1. Introduction 

In kinetic theory one often considers collisionless ensembles of classical particles 
which interact only by fields which they create collectively. This situation is 
commonly referred to as the mean field limit of a many particle system. More 
precisely such ensembles are described by a phase-space distribution function / : 
M x R 3 x M 3 — > Rj, where J D f(t, x, v)dx dv gives the number of particles which at 
time t have their position x and velocity v in the region D of phase-space K 6 . The 
Vlasov equation expresses the fact that / is constant along particle paths — which 
is a direct consequence of the absence of collisions — and reads 

d t f(t, x, v) + d x f(t, x,v)v + d v f(t, x, v) F(t, x) = 0, 

where F is some force, which will emerge via some field equation with its source 
given by the spatial particle density p(t, x) := f R3 f{t, x, v) dv. In the case of non- 
relativistic gravitational or electrostatic fields the corresponding system of partial 
differential equations is the Vlasov-Poisson system, in the case of relativistic elec- 
trodynamics it is the Vlasov-Maxwell system and in the case of general relativistic 
gravity the Vlasov-Einstein system. Such systems have been studied extensively in 
the literature; for an overview see [5l 116] . 

The system which is best understood is the Vlasov-Poisson system where p acts 
as a source term for the Poisson equation. For this system global-in-time classical 
solutions for general (compactly supported) initial data have been established in 
[T71 [P21 125] ; for a review see [TH QI5] . The existence-theory for the other systems 
mentioned above is not equally well understood. For the relativistic Vlasov-Maxwell 
system global classical solutions are so far only known in special cases (e.g. [7j 
Hn|)j while global weak solutions for general data were obtained in 0]. An 
investigation of the Vlasov-Einstein system was initiated in [ST] . 
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It is also interesting to note that the Vlasov-Poisson system possesses as (formal) 
singular limit cases the Euler-Poisson system with pressure zero and the classi- 
cal n-body problem, for both of which in general no global-in-time solutions ex- 
ist. More precisely, in the first case if one considers a phase-space density func- 
tion which is concentrated in u-space, i.e., f{t,x,v) — p(t,x)5(v — w(t,x)), where 
5 denotes the Dirac <5-function and w is a velocity field, then / formally solves 
the Vlasov-Poisson system iff (/, w) solves the pressure-less Euler-Poisson system. 
Similarly a density function / concentrated in position and velocity space, i.e., 
fit, x, v) = S(x — Xk(t))S(v — Vk{t)), formally solves the Vlasov-Poisson sys- 

tem iff Xk, Vk solve the n-body problem. One main problem here is — of course — the 
use of distributions in the context of nonlinear equations. Only few rigorous results 
relating (approximating sequences of) such concentrated solutions of the Vlasov- 
Poisson system to the solutions of the respective limit systems have been achieved; 
see [22] for the first case and [13] for the second. The main interest of course would 
be to use the far better existence-theory in the case of the Vlasov-Poisson system 
to learn something about the solutions of the related systems, e.g., in the context 
of shell crossing singularities in the case of the pressure-less Euler-Poisson system. 
In this work we propose the use of algebras of generalized functions (in the sense of 
J.F. Colombeau [2J[3]) to study these singular limits of the Vlasov-Poisson system. 
As a first step we prove an existence and uniqueness result for singular solutions, i.e., 
solutions concentrated either in position-space or in momentum-space (or both), to 
the spherically symmetric Vlasov-Poisson system in a suitable algebra of general- 
ized functions, where the latter provides us with a consistent framework for treating 
singular i.e., distributional solutions of nonlinear PDEs. The fundamental strategy 
of solving PDEs with singular initial data in the setting of algebras of generalized 
functions is regularization of singularities by convolution with a mollifier depend- 
ing on a regularization parameter e and first solving the equation for fixed e using 
existence theory in the smooth setting. Proving existence and uniqueness of gen- 
eralized solutions then amounts to deriving asymptotic estimates with respect to 
the regularization parameter. This process may alternatively be seen as uncovering 
new asymptotic stability results of smooth solutions to the system under pertur- 
bations of the initial data, which in our view is of independent interest. For a 
general discussion of applications of Colombeau theory to PDEs see [15]. Recent 
investigations into linear PDEs in this framework can be found in [9] HOj . 
We organize our presentation in the following way. In section [2] we collect some 
well-known facts on the spherically symmetric Vlasov-Poisson system which will 
be used later on and recall the basic definitions of generalized function algebras in 
the sense of J.F. Colombeau. Our main results are stated and proved in section [3] 
Finally we collect some facts on solutions of the Poisson equation in this setting of 
generalized functions in an appendix. 

Although our notation is mostly standard or self-explaining we explicitly mention 
the following conventions: For a function h = h(t, x, v) or h = h(t, x) we denote for 
given t by h(t) the corresponding function of the remaining variables. By || . || p we 
denote the usual L p -norm for p G [1, oo]. The index c in function spaces refers to 
compactly supported functions. Constants denoted by C may change their value 
from line to line but never depend on e. 
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2. Preliminaries 

We start by collecting some preliminaries (for a comprehensive presentation includ- 
ing full proofs see [TH]) from the existence-theory of the (spherically symmetric) 
Vlasov-Poisson system, which from now on we shall abbreviate by (VP) 

dtf + vd x f -d x ud v f = (1) 
Am = 4n-fp (2) 



P = fdv (3) 
where 7 = ±1. We suppose the following initial rcsp. boundary conditions 

f(0,x,v) = f(x,v)>0eC^(R 6 ) (4) 
lim u(t,x) = 0. (5) 

\x\ — >oo 

We shall often combine position and velocity into a single variable z = (x,v) and de- 
note by Z(s) = Z(s, t, z) — (X(s, f, z), V(s, f, z)) the solutions of the characteristic 
system of ([1]), 

X(s) = V(s) 
V(s) = ~d x u{s,X(s)) 
with initial condition Z(t,t,z) = z. The solution of the Vlasov equation is then 

o 

given by f(t, z) = f(Z(0,t,z)), hence all i p -norms (1 < p < 00) of / are constant 
in time, as is the Z^-norm of p, i.e., the mass, which will be denoted by M. Clearly 
f(t) is compactly supported and we denote its velocity support by P(t) :— sup{|i;| : 
(x,v) £ supp(/(f))}. 

We shall call a function j : I 6 = M 3 x I 3 ^ 1 spherically symmetric if for all 
AeSO(3) 

g{Ax,Av) = g(x,v). (6) 

o 

It is well known that in case the initial value / of (VP) is spherically symmetric the 
respective solution /(f) will also have this property. Moreover, the spatial density 
p(t) will be spherically symmetric (in the usual sense on K 3 — we shall denote it 
hence by p(t,r), where r = \x\) and the Poisson equation simplifies to 

Au(f,r) = ^r 2 n'(t,r)j =A^p{t 1 r). 

By a slight abuse of notation, in what follows we will use u(t,x) and it(t, r) inter- 
changeably. In addition to the usual key-estimate on the solution of the Poisson 
equation with compactly supported source term p(t), i.e., 

\\dMmo <c\\ P (t)\\i /3 \\p{t)\\y 3 (?) 

in the present setup we also obtain the estimates 

M 

\d x u(t,r)\ < -y and (8) 

||^ +ei+ej «Wlloo < C||«*)IU («£N 3 u = l,2,3). (9) 
Combining equations and (JS]) one obtains for all r > 

\d x u(t,r)\ <Cmin(i P(t) 2 ). (10) 
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Note that the latter estimate together with the fact that from £ £ C 2 ([0, t]), g £ 
Z X (R) and ||(s)| < £f(f(s)) VO < s < f it follows that \i(t) - f(0)| < 2\/2||c;||J /2 
yields boundedness of P(t) (cf. [IS], proof of Thm. 1.4). Hence global existence of 
solutions follows by the standard continuation criterion. 

We now turn to algebras of generalized functions in the sense of J.F. Colombeau 
[21 [3] . These are differential algebras containing the vector space of distributions 
V as a subspace and C°° as a subalgebra while displaying maximal consistency 
with respect to classical analysis according to L. Schwartz' impossibility result [24] . 
The main ingredient of the construction is regularization of distributions by nets of 
smooth functions and asymptotic estimates in terms of the regularization parameter 
e £ (0, 1] =: /, which in our case will be L°°-estimates global in z on compact time 
intervals. We shall work within the so-called special version of the theory and use [H] 
as our main reference. Colombeau algebras are defined as quotients of the spaces 
of moderate modulo negligible nets (u e ) s in some basic space £. In the present 
case we use £ = C°°(Rq x W 1 ) 1 and the following estimates for moderateness and 
negligibility (where O denotes the Landau symbol). 

£ S M {R+ x R") := {(u E ) E ££: VK CC Rq" Vch £ N^ l+1 3N £ N : 

sup \d a u £ (t,z)\ =0{e- N ) (ase^O)}, (11) 

(t.z)eKxRn 

JVg(R^ x K") := {(u e ) e ££: VK CCR+ Va£ N' l+1 Vm £ N : 

sup \d a u £ (t,z)\ =C>(e m ) (ase-+0)}, (12) 

Gg(R+ x R") := 4 f (R+ x M")/AAg(M+ x R n ). 

The index g in the above definitions signifies the global estimates with respect 
to x, v (contrary to the local estimates w.r.t. t). Generalized functions shall be 
denoted by u = [(u e ) e ], meaning that u is the equivalence class of the net (u e ) e . 
In the following section we shall prove existence and uniqueness results for the 
spherically symmetric (VP)-system in this setting. By a solution of a differential 
equation in a Colombeau algebra Qg we mean an element [(it E ) e ] of the algebra 
such that each u e solves the equations up to an element of the ideal A/g . Roughly, 
establishing solvability of a PDE in this setting therefore amounts to obtaining a 
classical solution u e for each e and proving moderateness of the resulting net (u £ ) e . 
Proving uniqueness amounts to showing that any two nets solving the equation up 
to an element of the ideal (with initial data differing by an element of the respective 
ideal) necessarily belong to the same equivalence class in Q g (cf., e.g., [15]). 
We note the following fact which will be essential in the uniqueness-part of the 
proof of our main result (and follows by an easy adaptation of the proof of [8] Thm. 
1.2.3]): u = [{u e ) e ] £ fffQEtf x R") is negligible iff 

VK CC Rj Vm £ N : sup \u £ (t, z)\ = 0(e m ) . 
We shall also need a suitable algebra of generalized functions containing the initial 

o 

data /. To this end we consider nets (u e ) E £ C°°(R n ) satisfying estimates of the 
form 

Va £ Nq 3N £ N (rcsp. Vm £ N) : sup \d a u £ (x)\ = 0{ £ - N ) (resp. 0{e m )). 
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We denote the respective spaces by £ 9 M , N g and Q g (cf. [I|fl4]). 
A function u in Q g (R n ) is called compactly supported if there exists a representative 
(it E ) £ of u and a compact set L containing the supports of all u £ . In this case we 
call the representative (u £ ) £ compactly supported. Note however, that since Q g is 
not a sheaf there is no well-defined notion of support for its elements (see Example 
below). 

The space T>' Lac (R n ) of bounded distributions (distributional derivatives of bounded 
functions) can be embedded into <? g (R™) by the map 

w i-> [(w * ip £ ) £ ] 

where ip is a rapidly decreasing function with unit integral and all higher order 
moments vanishing, and f £ (x) — e^ n ip(x/s). This embedding commutes with 
partial derivatives. Analogously, C°° (Rj , T)' LOO (R™)) can be embedded into Q g (R.q x 
R") via convolution. 

3. Generalized solutions of the Vlasov-Poisson system 

In this section we will state and prove our main results, providing existence and 
uniqueness of generalized solutions of the spherically symmetric Vlasov-Poisson 
system. We begin with a discussion of the relevant symmetry properties. 
We will call a generalized function g £ (? S (R 6 = R 3 x R 3 ) spherically symmetric if it 
possesses a representative (g e ) e that is spherically symmetric in the sense of ([6|) for 
all e. Likewise we call a function g £ C7 g (R 3 ) spherically symmetric if it possesses a 
representative (g e ) e that for fixed e is spherically symmetric in the usual sense. 
The following definition singles out classes of scales which can be used to measure 
the 'maximal degree of divergence' admissible in the initial data of (VP) to allow 
for unique solvability in the Colombeau algebra: 

Definition 3.1. Let p > 0. 

(i) By Sp 1 ' we denote the space of all scales a : I 
e — > and 

o-(s)- 1 = Odlog^l 1 ^) (e- 

(ii) By Sp 2 -* we denote the space of all scales a : I 
e — > and 

VC* > exp(C*a(e)- p ) = 0{\ log(e)|) (e -»■ 0) . 

Then Sp^ C Eg for p > q and i = 1,2. Note that using any scaling a sat- 
isfying <j(e) — > for £ — » 0, a 5-source can be viewed as the element [(<p a ( e ) ) £ ] 
of the Colombeau algebra. Since obviously ip„ E — > <5 in T>' as e — > 0, any such 
delta net is associated to the standard image [(y e ) e ] of the Dirac measure, hence 
macroscopically indistinguishable from it (cf. 3, 15, 8J for discussions of the concept 
of association and its effects on nonlinear modelling of singularities). After these 
preparations we may state our existence result. 

Theorem 3.2. (Existence of generalized solutions) 

o 

Let f £ Q g (R ) with a spherically symmetric, non-negative and compactly supported 

o 

representative (f e ) £ satisfying 

o 

(i) ||/e||i = M (the mass) for all e, and 



-> / satisfying a(e) — * for 
0). 

-> / satisfying a(e) — > for 
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( 1) ° 

(ii) there exists some a £ Yy 2 such that ||/ e ||oc < 



C 



(7(e) 

Then there exists a solution (f,u) of (VP) in ^(Rq x r6 ) x £g( K o~ X R 3 ) with 

o 

f(0,x,v) = f{x,v) and u vanishing at infinity (in the sense of Definition \A.5\) . 
Uniqueness of generalized solutions needs stronger assumptions on the data. We 

(2) 

present two results; the first one requires a E), -scale. 

Theorem 3.3. (Uniqueness of generalized solutions in Q~ g ) 

Let the assumptions of Theorem \3.2\ be satisfied but strengthen (ii) to 



(2) ° C 

(ii') there exists some a £ Yr 2 such that ||/ e ||oo < 



o{e) 

Then (/, it) given in Theorem ] 3. OH is the unique solution of (VP) with f(0,x,v) = 

o 

f(x,v), velocity support of f bounded as in Lemma Iff. 51 (ii) and u vanishing at 
infinity (in the sense of Definition \A.'E\) . 

If we adjust the algebra to the symmetry of our problem we can do without a 
Ep -scale. More precisely we change the basic space S in the definitions of £ a Ml 
N a and g 5 , respectively to E° := {(u e ) e G S = C°°(M^ x W 1 ) 1 ] Vi G M+ Ve £ 
I : u £ (t) is spherically symmetric }, where in case n — 6 spherical symmetry is 
to be understood in the sense of ([6]) and in case n — 3 in the usual sense. We 
denote the resulting algebra by . Likewise in case of the algebra Q g we take nets 
(u £ ) £ £ C°°(W 1 ) 1 such that u e is spherically symmetric, again in the respective 
senses for n = 6 and n — 3. The resulting algebra is denoted by Q°. Now we may 
state. 

Theorem 3.4. (Uniqueness of generalized solutions in Q?) 

Let f £ Q°(M. ) with a compactly supported and non-negative representative (f e ) e 
satisfying (i) in Theorem \3.2\ and 

( l) ° C 

(ii") there exists some a £ E 10 such that ||9"/ e ||cx) < Y+\a\ ^ 0T ' a < — ^' 

Then there exists a unique solution (f,u) of (VP) in G~(R$ x R 6 ) x G~{M£ x M 3 ) 

o 

with f(0,x,v) — f(x,v), velocity support of f bounded as in Lemma \3.5\ (ii) and u 
vanishing at infinity. 

Note that the assumptions in the above theorems in particular allow to model 
concentrated data which lead to the singular limits of the Vlasov-Poisson system 
described in the introduction. 

To prepare the proof of Theorem 13.21 first note that for fixed e the classical theory 
provides us with a unique solution (/ £ ,u £ ) in C°°(]Rj,R 6 ) x C°°(Rq x M 3 ) with 

o 

initial data / e (0, z) = f e (z) and u e (t) — > as |x| — » oo. Moreover, the solution will 
inherit the symmetry property of the data, that is f e (t), Pe{t) as well as u e (t) will 
be spherically symmetric. 

To prove the existence of generalized solutions we have to verify the moderateness 
estimates in (TTTj) . We split this task into two Lemmas collecting the necessary 
estimates. 

Lemma 3.5. (Zero order estimates) There exists C > such that for e sufficiently 
small we have for all t £ R + 
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(i) ||/ e (t)||oo<^, ||/ E (i)||x = ||p £ (t)||l=M 
(ii) Pe(t) < -^r- 

c 

(Hi) |K(t)||oo, ||0xt«e(f)||<x> < — -7T 

cr(£)3 

(iv) ||PeW||oo<-^. 

Moreover, for any T > and e sufficiently small 
(v) sup te[0 . T] ||Z £ (i)||oo < — ^r. 

cr(e)3 

Proof, (i) The L°°-estimate follows easily since for e fixed by the smooth result (see 

o 

Section 2) we have f e (t, z) — f £ (Z £ (0, s, z)). The L -estimates are immediate from 

assumption (i) in the theorem. 

(ii) We conclude from equation ([7]) and (i) 

R« £ (*)||oc < C\\p E (tM\\ Ps (t)\\2c < C^r. (13) 

We set 

, > • f 1 ( P e{t)\ 2 

g e {t,r) := mm — , I — — r 

Note that g £ (s,r) < g £ (t,r) for s < t since P e is monotonically increasing. Then 
combining \d x u £ (t,r)\ < M/r 2 with the above estimate we obtain from the charac- 
teristic equation 

\Xi(s)\ = \d Xi u E (*,Ms))\ < Cg e (t, \Xl(s)\) 

for s < t and 1 < i < 3. Therefore by the standard argument mentioned below 
equation (fTU|) we obtain 

|^(t)-^(0)|<2C>/2||<fe(t)||? 
and are left with calculating the L 1 -norm of g £ (t). We have 

OO OO -fe(t) OC 

J \g e (t,r)\dr = 2 J \g e (t,r)\dr < 2C J (^^dr + 2C J ±dr 

-oo g(e) i/3 

= 2^ + 2^<C F ^ 



Thus we obtain 



a(e)3 cr (£) 3 cr ( e ) 3 

\X*(t)-X*(0)\<C-±L- (14) 
ct(£) e 



and hence from the definition of P £ 

P e (t)<P+CP c (t)*ir(e)-i, 

o o 

where P bounds the diameter of the support of f £ . This in turn implies that P £ (t) 
is bounded independent of t for e fixed, and that P £ {t) < c T , which together 

<r(e)3 

with (fT4"|) gives (ii) and (v). 
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We insert (ii) into (fT3]) to prove (iii), i.e., 

C 



\\d x u £ {t)\\ c 



cr(e)3 

The estimate on u e (t) now follows easily by integration (taking into account that 
\u E (t, x)\ — 0(l/|x|)), while for (iv) we note 

C 



\p e (t)\\°o < cwfWMty < 



□ 



Lemma 3.6. (Higher order x,v- estimates) For all a £ N(j, all f3 € Nq and all 
T > there exists C > such that for e sufficiently small and all t £ [0, T] we have 

(i) ||^AW||oo< e CCT(Er2 

(ii) ||^p e («)|U < e c ^~ 2 

(iii) \\d?Z e (t)\\oo<e c °^- 2 

(iv) \\di +e ^u s (t)\\oo<e c ^ Vi,j. 

Note that compared with the zeroth order estimates we have to use an exponential 
term in a to bound the respective expressions necessitating the use of the scale 
a in condition [2] in Theorem 13.21 However, this term, i.e., exp(cr(e)~ 2 ) suffices to 
bound derivatives of any order. In particular, higher order derivatives do not lead 
to higher order exponential terms which would cause our approach to fail. 

Proof. We prove the Lemma by induction on \a\ and \/3\. 

In the case \a\ = \(3\ = we have shown even stronger estimates on f E (t), Z e (t) 
and p £ (t) already in Lemma l3.5l The only remaining estimate is the one on d 2 u E (t) 
which follows from ||0X(*)ll°o < C\\ Pe {t)\\™ < Ca{ £ y 2 . 

To carry out the inductive step we assume the Lemma holds for |a|, \P\ < n. We 
have to infer the respective estimates for \a\ = \/3\ = n + 1. We define 

:=a?x e (*), a«v s ( s ). 

Using the characteristic system we obtain (for suitably chosen i) 
£ a \s) = ±- s d«X e {s) = d«V e {s) = 4<*\ S ) 

= -8%~ e * (plu e (a, X e (s, t, z))d Zi X e (s, t, z)) 
= -d 2 x u E (s,X E (s))^(s) 

E ei )d2(d 2 x u E ( S ,X E (s)))dr~<Xs( S ) 

0<7<ct-ei \ ' / 

The last expression is a sum of products of terms of the form 

d x u E (s,X E (s)) with |i| < n + 2 

for which we have ||0*u e (*)||oo < C\\d s x ' p E (t)\\ OQ < exp(Ccr(e)- 2 ) by © and the 
induction hypothesis since \S'\ <n, and 

{d v z X E y' {d^X E y' with max(M, |w|) < n 
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which by induction hypothesis is also bounded by exp(Ccr(£)~ 2 ) on compact time 
intervals. 

So we find using © for |a| = and Lemma f3. 51 (iv) 

\4 a) (s)\ < \dZu £ (s,X E (s))\\^( S )\+e c ^)- 2 < Ca(e)- 2 \&\s)\ + e c ^~\ 
Hence summing up we obtain 

\4 a) (s)\ + \i { °Hs)\ < + Ca(e)- 2 {\4 a) (s)\ + \&H*)\), 

which by Gronwall's Lemma gives 

\4<*\s)\, \4 a) (s)\<e c ^r\ (15) 

i.e., \d?Z e (s)\ < exp(Ca{e)- 2 ) on [0,T] for all \a\ = n + 1, which is (iii). 
From here we obtain 

||a?/ e (t)||oo < Ce c °V~ 2 for all \a\ =n+l 

O j. o 

since d"(f e (t, Z s (0, t, z)j) is a sum of products of certain d z f £ (Z £ (0,t,z)) with 
products of powers of derivatives of Z £ (0, t, z), and we can use (iii) and the moder- 

o 

ateness of f £ . Thereby we have also shown (i). 

Item (ii) is now obvious using Lemma 13.51 (ii) . Finally, to prove (iv) we combine 
(ii) with (O. 

□ 

Proof Theorem\3M 

By the estimates of Lemma 13.51 and Lemma 13.61 above we obtain the necessary 
bounds on sup {t z)eKxRe \d^f E (t,z)\ and sup (t :!;)eK><R 3 \d^U e (t,z)\, where K is a 
compact subset of [0,oo). 

To obtain the estimates on dtf e we plug the estimates established so far into the 
Vlasov equation (using the bounded velocity support of f E (tj). From here the 
estimate on dtp £ and hence on dtu £ follows. Now differentiating the Vlasov equation 
we obtain the estimates on terms of the form dtd z f £ . Higher order dt- and mixed 
(t, z)-estimates of f £ are obtained by successively differentiating Vlasov's equation 
and in turn imply the respective estimates on u £ . 

Moreover u(t) = [(u £ (i)) e ] is vanishing at infinity in the sense of Definition IA.5I 
since the support of Au £ (t) is bounded by C + tP £ (t) and P £ (t) < Ccr(e) -1 / 3 by 
Lemma 13.51 (ii). 

This proves existence of solutions in Gg(R + x K 6 ) x Qg(R+ x R 3 ) with f(0,x,v) = 

o 

f(x,v) and u(t) vanishing at infinity. □ 
Proof Theorem^ 

We have to prove uniqueness of the solution obtained above. So assume (/ = 
l(f £ ) £ ],u — {(u £ ) £ ]) is a solution as constructed above and let (/ = [(f £ ) £ ],u — 
[(u £ ) £ ]) be another solution of the (VP) system with the same initial data (i.e., 
/(0) = /), u vanishing at infinity in the sense of Definition [A3] (with distinguished 
representative u £ ), and e-wise bounded (by P s {t), satisfying (ii) of Lemma |3.5|) 
velocity support of f £ (t). Proving uniqueness in our setting amounts to establishing 
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that the differences f £ — f s and u £ — u £ lie in the respective ideals. We have 
dtfe + vd x f £ ~ d x ii £ d v f £ = n e 

Au £ = 47T7 / f £ dv + n £ (16) 



/ 8 (0) = f e + n e =:f°, 

where (n £ ) £ denotes a "generic" (analogous to the "generic" constant C) element 
of the ideal which may denote different negligible quantities in each equation. De- 
noting by Z £ the characteristics of the above inhomogcneous Vlasov equation we 
obtain (cf. e.g., [11], appendix A): 

f £ (t,z) = f°(Z £ (0,t,z)) + f n e (s,Z e (s,t,z))ds 

Jo 

= f°(Z e (0,t,z))+n £ (t,z) = f £ (Z £ (0,t,z))+n e (t,z). (17) 

Consequently we may estimate the difference in the distribution functions using 
Lemma [531 (i) 

\f s (t,z)-f s (t,z)\ < \\dJ e \\ oo \Z s (0,t,z)-Z e (0,t,z)\ + \n e (t,z)\ 

< e c ^ 2 \Z £ (0,t,z)-Z £ (0,t,z)\ + \n £ (t,z)\. (18) 

For the characteristics we obtain (for < s < t) 

t 

\X e (a) - X e (a)\ < [ \V £ (s')-V £ (s')\ds> 



i 

\V £ (s)-V £ (s)\ < f\d x u £ (s',X £ (s'))-d x u e (s',X £ (s'))\ds' 

S 

t 

< J (\d x u £ (s',X £ ( S ')) - d x u £ (s',X £ (s'))\ 

S 

+ \d x u £ ( S ',X £ (s'))-d x ii £ ( S ',X £ ( S '))\) ds' 
t 

< sup feOOlU f\X e (s')-X e (s')\ds' 

s 

\d x u £ (s', X £ (s')) - d x u £ (s', X £ (s'))\ds. (19) 

Now we turn to the perturbed Poisson equation for u £ , i.e., 

Au e (t, x) = 47r 7/ 5 £ (i, x) + n £ {t, x). (20) 

Since u is strongly vanishing at infinity the right hand side in the above equation 
has its support in B £ -n(0) for some N > 0. Furthermore by assumption f £ (t) 
has its w-support contained in some Bp n\(0) and hence its a;-support bounded 

O £ ~ O 

by R+ L P e (s)ds, where R bounds the x-support of / e (0). This implies that the 
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support of f e (t) is bounded by some with Q £ (t) < C(t)P £ (t) where C(t) 

depends linearly on time. 

As a consequence n £ (t) in equation (|20[) above has its support also contained in 
some B £ -n(0). Therefore we may define n £ (t,x) := J n £ (t,y)/\x — y\ dy, which is 
clearly in the ideal and finally we have found a representative u £ :— u £ — n £ of [(£i e ) e ] 
that satisfies the non-perturbed Poisson equation with source p £l i.e., Au £ — An^p £ . 
This in turn implies A(u £ — u £ ) — 4:irj(p £ — p £ ) and using JJ]) we write 

\\d x (u e - u £ ){t)\\ 

<C\\ [ (f e (t,,v)-fS,,v))dv\\l /3 || / (/ £ (t,., w )-/ £ (<,., w ))d v ||^ 3 . 



On estimating the I/ 1 -norm above we use Q £ {t) := max(Q E (t), Q s (t)), where Q £ (t) 
denotes the respective bound on the support of f £ (t) and write 

|| f(f E (t,.,v) -f e (t,.,v))dy\ \x < [ \f £ (t,z)-f £ (t,z)\dz 



< CQ e (t) b \\f c (t)-f e (t)\\co 

< ^\\f £ (t)-f £ (t)\u 

for t £ [0, T], where we have used Lemma T3. 51 (ii) in the last step. So we find 

\\d x (u £ - U e )(t)||oo < ~°\\f e {t) ~ fe(t)\U (21) 

Now applying Gronwall's lemma to (fT§|) we obtain for all q 

2 / ' 

sup ||0 x « e (s )||oo r 

\Z £ {s) - Z £ {s)\ < Ce>^ \d x u E (s',X s (s'))-d x u e (s',X e (s'))\ds 



t t 
< e c ^ 2 I \\d x u £ {s') - fl^aOHoods' + / |M«')lloo<k' 



e c 



t 

° {er \^W^ J mS ' ] " ^')\\oods> + e«), (22) 



where we have again used Lemma 13.51 and (|2"T|) above. Finally we combine (fT5|) 
with ([22]) and use Gronwall's lemma for the second time to obtain for all q 

sup ||/ £ (t) - Ut)\\oo < exp (a-^e c ^ 2 ) e\ 
te[o,T] 

and due to our assumptions on the scale we see that the difference of the distribution 
functions is in the ideal. From here the respective estimates on the difference of 
the spatial densities and on ||M e (t) — M e (i)l|oo follow easily. □ 

Proof of Theorem\^\ 

As for existence just observe that (ii") implies (ii) and that by classical theory the 
solution inherits the respective symmetry properties of the data. 
To prove uniqueness we assume that (/ = [(f £ ) £ ],u — [{u £ ) £ }) G £/|(Rq x R 6 ) x 
£|(Rq x R 3 ) is a solution as constructed above and (/ = [(f £ ) £ ] 7 u — [(u £ ) £ ]) 
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is another such solution with the same initial data, u vanishing at infinity (with 
distinguished representative (u E ) s ) and the velocity support of f e (t) bounded by 
P E (t) satisfying (ii) of Lemma 13.51 We follow the proof of Theorem 13.31 up to 
estimate (1191) but now using spherical symmetry we provide a stronger estimate on 



\d x u E (s',X E {s')) - d x u E (s' 7 X E (s'))\ds' 



< 4tt 



Xs(s') 



fl(s') 



fe(s') 



l 2 \ Pe {s',l)- p E (s',l)\dlds', (23) 



where r E denotes the modulus of X E . Note that this formula does not hold unless 
we use the algebra £/| since in general p E (t) will not be spherically symmetric due to 
the non-symmetric perturbations in (|16p . Estimating the difference of the spatial 
densities we find using P E (t) — max(P e (i), P E (t)) as well as (|17j) 

\p E {t,r)-p E (t,r)\ < t \f E (t,r,v)-f E {t,r,v)\dv 



< 



c 

a(e) 



\f e (Z e (0, t, z)) - f e (Z E (0 7 1, z)) I dv + \n e (t, r)\ 



lia.MU ||^(0, t, .) - Z s (0,t, .)|U + \n £ (t,r)\. 



Inserting this into (|23p we obtain using Lemma |3. 5 



u e {s',X e (s'))-d x u e {s',X E (s'))\ds 



<4"T H^Alloo sup \X e (s')\ I {}\Z e (0, S ',.)-Z e (0,s',. 
cr(e) s <s'<t 



+ sup \n e (s' , r)|)ds' 
<C*(e)-^ I \\Z s (0,s',.)-Z e (0,s',.)\\ oo ds' + Cs q 



for t G [0, T) and all q. Now combining this with we obtain 



|Z e (0,a ) .)-lr e (0,s ) .)||oo < Ca(e)~- / ||Z E (0, s', .) - Z s (0, s', .)||<x»*»' + Ce q 



which by Gronwall's lemma gives 



sup \\Z E (0,t,.)- Z E (0,t,. 
te[o,T] 
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Hence by our assumption on the scale the difference of the characteristics is negli- 
gible. By (fTT|) this immediately implies that [(/ e )J = [(/ £ )e] and so the same holds 
true for the spatial density as well as for the potential. □ 

Appendix A. Uniqueness for generalized solutions of the Poisson 

EQUATION 

In this appendix we collect some facts on the Poisson equation within the frame- 
work of nonlinear generalized functions. We focus on the question of uniqueness, 
presenting a solution concept providing the existence of unique generalized solu- 
tions subject to a boundary condition generalizing the classical condition u — > 
(| a; | — > oo). Throughout this appendix we assume that n > 3 and write the Poisson 
equation as Au = p. Also, we denote the fundamental solution of the Laplace 
equation by C n /\x\ n ~ 2 . 

In addition to the algebra Q g (R n ) used in our main results we also treat the case 
of the standard (special) Colombeau algebra <7 s (f2) (with C R") which is defined 
using estimates on compact subsets of J7, i.e., 

Va G N' l VK CC n 3N G N (resp. Vm G N) : 
sup \d a u £ (x)\ = 0(e- N ) (resp. 0(e m )). 

xEK 

We begin with some preliminaries. Let u G Q s (fl). Then u has compact support 
(that is: 3K CC fl : u\q\k = 0) if and only if there exists a representative (u s ) e of 
u and L CC £1 such that supp(w e ) C L for all e > 0. In this case we say that (u e ) e 
itself has compact support. 

Indeed for any compactly supported u we may choose a cut off function \ G T>(fl) 
such that y = 1 on a neighborhood of the support of u. Then for any representative 
{u e ) £ of u we construct a new representative (yu £ ) e which vanishes outside the 
support of x- 

However, in general L will properly contain the support of u in its interior. Indeed 
let u = l(5) (with i denoting the embedding of distributions into the algebra of gen- 
eralized functions) then there clearly exist representatives that vanish outside any 
compact neighborhood of the origin. On the other hand there is no representative 
which vanishes outside the support of u. 

Next we note that any generalized function u G £ s (]R n ) has a representative (u e ) e 
which vanishes at infinity, i.e., u e (x) — ► (|x| — » oo) Ve. Indeed take any represen- 
tative of u and multiply it with an e-dependent cut-off function \e which is equal 
to unity inside a ball of radius l/(2e) and vanishes outside a ball of radius 1/e. 
Moreover we have the following warning example of non-uniqueness of generalized 
solutions to the Laplace equation. 

Example A.l. We consider Au = in C/ s (K") . Clearly u = is a solution. On 
the other hand we construct a solution u as follows: Set u e = \e with Xe & s above. 
Then ii E vanishes at infinity, [(u e )] — 1 a> n d Au = 0. 

However, there does not exist a representative (u £ ) e of u such that supp(AM e ) is 
contained in some ball of radius R for all e. Indeed suppose to the contrary that 
Aii £ = n £ G M S (W 1 ) with supp(n e ) G Br(0) for all e. Then by classical uniqueness 
we have that u e {x) — C n J ^J^^-'i dy and hence (u £ ) s is in the ideal which is not 
possible. 
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This observation motivates the following definition securing uniqueness of solutions 
to the Poisson equation. 

Definition A. 2. Let p 6 C* S (R™) be compactly supported. We call u € Q s {W l ) a 
solution of the Poisson equation vanishing at infinity if Au — p and if there exists 
a representative (u £ ) £ of u that satisfies 

(i) Ve > : lim u £ (x) = 0, and 

x — >oo 

(ii) (Au e ) £ is compactly supported. 
We may now state the following result. 

Theorem A. 3. Let p G £/ s (R n ) be compactly supported. Then there exists one and 
only one solution of the Poisson equation 

Au = p 

vanishing at infinity. 

Note that the assumptions on u in Definition I A. 2 1 arc not redundant. Indeed the 
compact support of p guarantees the existence of a representative {u £ ) £ satisfying 
property (ii) and there also exists a representative (u e ) £ of u which vanishes at 
infinity. However, in general u £ ^= u £ and uniqueness may fail as is explicitly 
demonstrated by the example above. 

Proof. Existence: By the above we may choose a compactly supported representa- 
tive (p £ ) £ of p and define 

u £ {x) :=C n [ P<i{V l_ 2 dy 
J \x-y\ n 1 

By the classical theory u £ satisfies both requirements stated in the theorem. 
Uniqueness: Let u, u be two solutions as above and choose representatives {u £ ) £ 
and (u e ) s satisfying (i) and (ii) in Definition IA.2I From the second property we 
conclude that A(u £ ~ u £ ) = n £ is compactly supported. By the first property 
we have u £ — u £ — > (|x| — > oo). Hence by the classical theory (u £ — u £ )(x) — 
C n J \J^y\-}-2 dy which obviously is in the ideal. □ 

We now turn to the "global" algebra Q g . The basic difference between Q g and Q is 
that due to the global estimates defining it, Q g is not a sheaf: 

Example A. 4. Define u £ € C°°((— m, m)) to be 1 for e > l/m and exp(— 1/e) 
fore < l/m. Choose a partition of unity (xm)meN subordinate to ((— m, m)) m£ N 
and set 

oo 

u e (x) := J2xm(x)4 m Hx). 

m—1 

Then (u £ ) £ £ ££f(R) \J\f g (M), so u — [(u e ) s ] provides an example of a nonzero 
element ofQ g (M) whose restriction to each (— m, m) is zero. 

Clearly in this setting Example I A . 1 1 does not work since here (u £ ) £ is not a repre- 
sentative of the function 1. This opens the possibility of relaxing condition (ii) in 
Definition IA. 21 which is necessary in the context of the (VP)-system since Au £ {t) as 
constructed in the proof of Theorem 13.21 will not be compactly supported. On the 
other hand we have proved P £ {t) < Ccr(e)^ 1 / 3 in Lemma T3. 5 1 (ii). This motivates 
the following definition which will provide us with the solution concept used in our 
main results. 
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Definition A. 5. Let p € C7 g (R n ) be compactly supported. We call u £ Q g (M. n ) a 
solution of the Poisson equation vanishing at infinity if Au = p and if there exists 
a representative (u £ ) £ of u that satisfies 

(i) Ve > : lim u e (x) — 0, and 

x — >oo 

(ii) supp(Aw e ) e C B £ -n(0) for some N > 0. 

Note that again conditions (i) and (ii) are not redundant. Indeed take u £ with 
u £ = 1 on _B e i/ e (0) and vanishing outside a ball of twice that radius. Then (i) 
clearly holds but Au £ ^ near |af| = e 1/,£ . The desired result in this framework is 

Theorem A. 6. Let p 6 Q g (M. n ) be compactly supported. Then there exists one and 
only one solution of the Poisson equation 

Au = p 

vanishing at infinity. 

Proof. Existence is proved as in Theorem IA. 31 

To prove uniqueness suppose we have two solutions u, u in Q g (R. n ) vanishing at 
infinity. Let (u e )e an d (u e ) e be representatives according to Definition IA.5I By 
condition (ii) we have A(u £ — u s ) = n £ with (n e ) e in the ideal and supp(7i e ) C 
B e -n(0) for some N. So 

— N 

I u £ - u £ \(x)<C n I , W[y }}_ 2 dy < Ce m [ rdr Vm, 

Jb c - n (o) \ x — y\ Jo 

hence is in the ideal. □ 
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